Critically balanced ion temperature gradient turbulence in fusion plasmas 
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Scaling laws for ion temperature gradient driven turbulence in magnetized toroidal plasmas are 
derived and compared with direct numerical simulations. Predicted dependences of turbulence fluc- 
tuation amplitudes, spatial scales, and resulting heat fluxes on temperature gradient and magnetic 
field line pitch are found to agree with numerical results in both the driving and inertial ranges. 
Evidence is provided to support the critical balance conjecture that parallel streaming and nonlinear 
perpendicular decorrelation times are comparable at all spatial scales, leading to a scaling relation- 
ship between parallel and perpendicular spatial scales. This indicates that even strongly magnetized 
plasma turbulence is intrinsically three-dimensional. 



Introduction. Microscale turbulence is a ubiquitous 
feature of the plasmas used for magnetic confinement 
fusion. It is driven by kinetic instabilities feeding pre- 
dominantly off a strong mean gradient in the ion tem- 
perature, and it is responsible for the majority of parti- 
cle and heat transport observed in experiment. As with 
neutral fluid and magnetohydrodynamic turbulence, ex- 
act analytical results for kinetic plasma turbulence are 
rare, and numerical simulations are costly. Phenomeno- 
logical scaling laws are thus useful for guiding simulation 
and providing gross predictions of plasma behavior in a 
multi-dimensional parameter space. 

Experimental, numerical, and analytical results (cf. 
[IHlj) have long been used to predict the dependence of 
turbulent fluxes on the mean plasma gradients and on the 
magnetic field configuration. However, scalings based on 
empirical observations provide limited physical insight, 
and the theoretical predictions, which are predominantly 
based on linear or quasilinear arguments, are not suffi- 
ciently detailed to be easily falsifiable. A more detailed 
examination of the properties of kinetic plasma turbu- 
lence has been conducted for scales smaller than the ion 
Larmor radius [SHl], but it is the ion temperature gra- 
dient (ITG) driven turbulence above the Larmor scale 
that is most relevant for heat transport in fusion devices 
(cf. [8]). Recent advances in plasma fluctuation measure- 
ments O [To have provided turbulence spectra in this 
scale range; direct numerical simulations have also calcu- 
lated spectra [TT] and studied energy injection, transfer, 
and dissipation [T2j [13] . 

In this Letter, we propose a phenomenological scal- 
ing theory of ITG turbulence. A number of simple, 
physically- motivated conjectures about the nature of this 
turbulence are formulated and applied to obtain fluctu- 
ation spectra from the driving scale to the ion Larmor 
scale. We then derive predictions for the dependence of 
heat flux on plasma current and ion temperature gra- 
dient. Numerical results are presented to support our 
predictions and justify our conjectures. 

Gyrokinetic turbulence. Plasma fluctuations in a 
strong mean magnetic field are anisotropic with respect 
to the mean field direction and have typical frequencies 



that are small compared to the ion Larmor frequency, 
Qi. Such fluctuations are correctly described by the 
gyrokinetic approximation It assumes rfti 1, 

~ Pi/L <C 1, and pi/i± ~ 1, where r is the fluc- 
tuation time, pi the ion Larmor radius, L the character- 
istic length scale of the mean dynamics, and i\\ and i± 
are the fluctuation length scales parallel and perpendic- 
ular to the mean field, respectively. Averaging over the 
fast Larmor gyration eliminates gyroangle dependence. 
For the remaining phase space coordinates, we choose 
(R, /i), where R is the position of the center of a par- 
ticle's Larmor orbit, E = is the kinetic energy of a 
particle, and /i = mv\/2B is the magnetic moment, with 
m particle mass, v particle speed, v±_ its perpendicular 
component, and B mean magnetic field strength. 

We restrict our attention to ion scale turbulence with 
Pi/^^ < I7 foi" which the guiding center position, R, and 
actual particle position, r, are approximately equal. The 
equation describing the evolution of electrostatic fluctu- 
ations in the absence of sonic flows is 

d ( 'Z's^^ \ ( \ 

dt ~ ~J^^^'' ) ^ ^ ^^''^ ' (1) 

^ ^^E {Fm^s ^ hs) = C[hs\, 

where 5fs = hs — {ZgCip /Ts)Fm,s describes the distribu- 
tion of particle positions and velocities for species 5, t 
is time, vm,s = (b/l^^) x (v|b • Vb + {v\/2)VB/B) is 

the magnetic drift velocity, v^; = {c/B)h x \/(p is the 
E X B drift velocity, b is the unit vector along the mean 
field, ZsC is the species charge, e is the proton charge, Tg 
is the temperature, Fm,s is a Maxwellian distribution of 
velocities, and (7 is a Fokker-Planck collision operator. 

The electrostatic potential cp is obtained by imposing 
quasineutrality, ZsCris = 0, which can be written 




where ris is species density. Assuming J d^vhi ~ '^th^*' 
where i denotes the main ion species and ^th its thermal 
speed, Eq. ^ gives hi/Fu^i ^ Zieip/Ti. 
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Critical Balance. Because the turbulence we are con- 
sidering is anisotropic, dimensional analysis alone is not 
sufficient to determine scalings of the fluctuation ampli- 
tudes with both and £_\_. To flx the ratio of to 
we make a conjecture known as critical balance fiE : The 
characteristic time associated with particle streaming and 
wave propagation along the mean field, ^y/'^th; is compa- 
rable to the nonlinear decorrelation time at each scale. 
This is motivated by the causality constraint: two points 
along the mean fleld can be correlated only if informa- 
tion can propagate between these points in the time it 
takes turbulence to decorrelate in the plane perpendic- 
ular to the fleld. It gives us a relation between parallel 
and perpendicular spatial scales: 



Vth 



(3) 



where Tni is the nonlinear decorrelation time, R is the 
major radius of the torus, = {eipi/T){R/ pi) and ipi = 
(P{t -\- i) — (p{r). The subscripts x and y refer to the 
coordinates in the plane perpendicular to the mean fleld, 
with X labeling surfaces of constant magnetic flux and y 
labeling fleld lines within a constant flux surface. 

Outer Scale. We deflne the outer scale as the scale 
for which the time associated with the linear drive is 
comparable to the nonlinear decorrelation time: 



'nl 



PiVth 



(4) 



where uj^ is the frequency associated with the ITG drive, 
Lt is the ITG scale length, and o labels outer scale quan- 
tities. The outer scale corresponds to the injection range, 
which contains the turbulence amplitude peak. We con- 
jecture that the characteristic parallel length at the outer 
scale, i^^, is the parallel system size. In toroidal plasmas, 
this is the distance along the mean fleld from the out- 
side to the inside of the torus. This distance, known as 
the connection length, is qR^ where the safety factor q 
measures the pitch of the fleld lines. Thus, 



q - qR. 



(5) 



Combining relations ([3|-([5| provides a prediction for the 
dependence of iy on q and R/Lt 120] 



Pi 



qR 



(6) 



Since qR/LT ^ 1 in fusion plasmas, the ratio pi/£y <C 
1. Thus, there is a range of scales between the outer 
scale and the Larmor scale, below which kinetic damping 
effects are expected to become signiflcant. 

The flnal piece of information necessary to determine 
scalings for the fluctuation amplitude at the outer scale, 
^o, is a relationship between and ly. We conjecture 
that the length scales in the perpendicular plane, ix and 
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FIG. 1: (a) Expectation value of kyPi versus qk,, where k, = 
R/Lt- (b) Normalized heat flux versus qk^ . Lines show the 
predicted scalings (|6| and 



iy, are comparable at all scales: £x ^ ~ Indeed, 
one might argue that is set nonlinearly at the outer 
scale through the shearing of radially extended eddies by 
zonal flow: ~ £y^{Sz¥^ni)i where Szy is the shearing 
rate at the outer scale due to zonal flow. For strong ITG 
turbulence, one expects SzF^ni ^ 1 [16 , so that ix ^ iy 
is satisfled. 

Taking ix ^ iy and using ^ and ([6| in (|4| gives 



pi Lt 



R 



(7) 



If we assume that the ion turbulent heat flux through 
volume V, Qi = V'^ J d^r J d^v {ve ■ Vx) {rrnv^ /2) Sfi, 
is dominated by the contribution from the outer scale, 
relations (|6| and ^ imply a scaling for Qi'. 



Q^ fR 

riiTiV th \Pi 



Pi 



R^ 



(8) 



The R/Lt scaling is only valid for sufficiently large R/Lt 
because our simple analysis ignores the flnite critical tem- 
perature gradient associated with the ITG instability. 

Inertial Range. We now consider the range of scales 
between i^ and pi and conjecture that it is an iner- 
tial range: There is no significant dissipation or driv- 
ing at scales between ij_ and pi. This conjecture is 
to be checked a posteriori. To determine the spec- 
trum in the inertial range, we identify the free energy 
W = Yjs I I d^^'^s^fs/^M^s as a nonlinear in- 
variant [5 and consider scale-by-scale energy balance. 
Because we are in an inertial range, the flux of free en- 
ergy, Wi/Tni^ must be independent of i_\_: 



1 We 



^M'^^^,'- constant. (9) 



^iTi Tfii \RJ R ixiy 

The dependence of ^£ on q and R/Lt is obtained by 
solving relation ([9| for matching to [see ([7|], using 
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FIG. 2: Electrostatic fluctuation spectra, E{ky), normalized 



to the predicted scaling ( 11 ) at the outer scale. We have nor- 
malized kyPi using the outer scale relationship (Ig]). Asterisks 
denote the Cyclone Base Case = 1.4, k = R/Lt = 6.9); 
circles and triangles denote simulations where q and k, vary 
from 2.8 to 7.0 and 10.0 to 17.5, respectively. The solid line 
gives the predicted inertial range scaling ( |11| . 



relation (pi), and assuming isotropy: 
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(10) 



Parseval's theorem relates ^£ to the Fourier coefficient 
giving a scaling for the ID fluctuation spectrum, 
Elky), defined so that / dky piE(ky) = V'^ J d^r 



E{ky)^kyP,\^J,f^q'/'{^ 



8/3 



-7/3 



(11) 



where k± = 27r/i±. Using (10) to evaluate r^^, we find 
^* ~ {^±/^±_y^^^ni^ J confirming that the drive is sub- 
dominant to the energy transfer in the inertial range [21 . 

Using relation (10) and applying critical balance (|3| 
gives the scaling of £11 with £±: 



qR 



Pi qR) 



4/3 



(12) 



Converting this into a scaling for we obtain the par- 
allel structure function ~ q{R / LtY / R) . 

Numerical results. We next compare our scaling re- 
lations with numerical results obtained using the gy- 
rokinetic code GS2 [17^. We restrict our attention to 
electrostatic fluctuations with perturbed electron den- 
sity drie/rie = e((/p — ^)/Te, with the overline denoting a 
flux surface average. The magnetic geometry is as in the 
widely-benchmarked Cyclone Base Case (CBC) [18]: un- 
shifted, circular magnetic-fiux surface, with r/i? = 0.18, 
s = dhiq/dhir = 0.8, and R/Ln = 2.2, where r is the 
minor radius of the flux surface, and is the density 
gradient scale length. The parameters q and hi = R/Lt 
were varied over several simulations to obtain numerical 




FIG. 3: Correlation function, C (|14|). The solid black line is 
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the critical balance scaling (12) 



scalings. These simulations employed a small amount of 
upwinding along the magnetic field and hyper-dissipation 
in k^ [19 , cutting off the fluctuation spectra at /c^Pi — 1. 

In Fig. [Ij we show the q and k dependences of the 
normalized ion heat flux, Qi, and of 



{kyPiY 



(13) 



which is a good measure of the outer scale, provided the 
spectrum is sufficiently steep. The simulations agree re- 
markably well with our predicted scalings (|6| and (|8| [22] . 

The one-dimensional spectrum, E{ky)^ is plotted 
vs. kypi in Fig. [2] for several q and values. At scales 
smaller than the outer scale, all spectra follow the same 



power law, which agrees with (11). 

Our predictions for the critical balance scaling ( 12 ) are 
tested in Figs. |3] and |4| The parallel correlation function. 



C(ky,M) 



_ Ek^ Hk,. ky, = 0)^*(fe„ ky, = AO) 



(14) 

is plotted in Fig. [3| for the simulation with q = 4.2 and 
= 6.9. Here is the poloidal angle of the torus so that 
^ i\\/qR^ with ^ = at the outermost point on the 
flux surface. Our prediction ( [12] ) is given by the black 
line, which fits the data. The (^i^-normalized correlation 
length, Ae{ky) = J d{M)C{ky,M), is plotted in Fig. [i] 
for multiple q and n values. At sufficiently large kypi^ the 
data follows the power law (12). Note that A6{ky) peaks 
at a value of approximately unity, in agreement with ([5| . 

Sub-Larmor scales. Using arguments similar to those 
presented above, scalings of E{k±) ~ kj_^^^^ and 
{k±Pi)c ^ Do^/^ were obtained for the sub-Larmor scales 
in [5 and verified numerically in [6 . Here the subscript c 
denotes the cutoff wavenumber, and the Borland number 



Do = {Tp 



[5]-[7] is the kinetic plasma turbulence 



analog of the Reynolds number, with r^. the nonlinear 
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FIG. 4: Scaling of poloidal correlation length with qnkyPi 
for several simulations with varying q and n = R/Lt- The 
dashed line indicates the critical balance scaling (12). 
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FIG. 5: Cartoon of the fluctuation spectrum from the outer 
scale, {k±pi)o [Eq. (|6|], to the dissipation scale, {k±pi)c [Eq. 
(15)]. Scalings for k±pi > 1 are taken from [5]. 



time at pi and ua the ion-ion collision frequency. Using 
relations (10) and ([3|, we find: 



{k^p^), - Do^/^ - q 



1/5 



R_ 



4/5 



Vth_ 

UiiR 
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(15) 



Combining the results of [5] with those given here pro- 
vides a complete picture of the spectra of ITG turbulence 
from the driving to dissipation scales, shown in Fig. [5] 

Discussion. The main results obtained in this Letter 
are: a scaling of heat flux ([8| and dissipation scale ([l5| 
with q and R/Lt] a power law scaling for the electrostatic 
fluctuation spectrum (11); and a relationship between 
parallel and perpendicular length scales ([T2I) 



The heat 

flux scaling, confirmed numerically (see Fig. [T]), has a 
stronger dependence on R/Lt than is usually assumed by 
reduced models for turbulent transport. While this has 
little effect for near-marginal plasma turbulence, it may 
be significant in the vicinity of steep gradient regions. 
The power law prediction for the turbulence spectrum 



(11), also confirmed numerically (see Fig. [2|, appears to 
be consistent with recent experimental fluctuation mea- 
surements [10 . Along with the predictions for the dis- 
sipation scale (15), the spectrum could be used to de- 



sign large eddy simulations for gyrokinetics [12 and to 
guide resolution choices for direct numerical simulations. 
The critical balance conjecture (|3| has proven to be ro- 
bustly satisfied for the ITG turbulence considered here 
(see Figs. [3] and [4|. This indicates that ITG turbulence 
is an inherently three-dimensional phenomenon. 
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